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Abstract. We investigate large time existence of solutions of the Navier-Stokes-Boussinesq 
equations with spatially almost periodic large data when the density stratification is suffi- 
ciently large. In 1996, Kimura and Herring [T5] examined numerical simulations to show a 
stabilizing effect due to the stratification. They observed scattered two-dimensional pancake- 
shaped vortex patches lying almost in the horizontal plane. Our result is a mathematical 
justification of the presence of such two-dimensional pancakes. To show the existence of 
solutions for large times, we use £ -norm of amplitudes. Existence for large times is then 
proven using techniques of fast singular oscillating limits and bootstrapping argument from 
a global-in-time unique solution of the system of limit equations. 

1. Introduction 

Large-scale fluids such as atmosphere and ocean are parts of geophysical fluids, and the 
Coriolis force due to the earth rotation plays a significant role in the large scale flows con- 
sidered in meteorology and geophysics. 

Mathematically, it was first investigated by Poincare [19]. Later on, the problem of strong 
Coriolis force was extensively studied. Babin, Mahalov and Nicolaenko (BMN) [TJ [2] stud- 
ied the incompressible rotating Navier-Stokes and Euler equations in the periodic case while 
Chemin, Desjardins, Gallagher and Grenier [8] analyzed the case of decaying data and more 
recently, the second author [5D] considered the almost periodic case. Gallagher in [TO] studied 
a more abstract parabolic system. We also refer to Paicu [TO] for anisotropic viscous fluids, 
Benameur, Ibrahim and Majdoub [5] for rotating Magneto- Hydro-Dynamic system and to 
Gallagher and Saint-Raymond [TT] for inhomogeheous rotating fluid equations. 
Moreover on the one hand, the case when fluids are governed by both a strong Coriolis force 
and vertical stratification effects was investigated by BMN in [3] in the periodic setting and 
Charve in [6] for decaying data. However, their studies do not cover the case when fluid 
equations are governed by the only effect of stratification. It is known that a strong Coriolis 
force has a stabilizing effect (see [1]). However, in BMN [H Section 9.2] the authors observed 
that for ideal fluids (i.e., with zero viscosity), the only effect of stratification leads to un- 
balanced dynamics. Moreover, the case of both strong Coriolis and stratification forces in 
the almost periodic setting seems to remain open. Finally, note that for the almost periodic 
case, energy type estimates cannot be used, and instead Fujita-Kato's approach has to be 

S. I. is partially supported by NSERC# 371637-2009 grant and a start up fund from University of Victoria. 
T. Y. is partially supported by PIMS Post-doc fellowship at the University of Victoria, and partially 
supported by NSERC# 371637-2009. 



2 



SLIM IBRAHIM AND TSUYOSHI YONEDA 



used. On the other hand, Kimura and Herring [16] examined numerical simulations to show 
a stabilizing effect due to the effect of stratification for viscous fluid. They observed scattered 
two-dimensional pancake-shaped vortex patches lying almost in the horizontal plane. Our 
result can be seen as a mathematical justification of the presence of such two-dimensional pan- 
cakes. More precisely, we study long-time solvability for Navier-Stokes-Boussinesq equation 
with stratification effects. The Navier-Stokes-Boussinesq equations with stratification effects 
are governed by the following equations. 

' d t u - uAu + (u • V)u + X7p = gpe 3 x G M 3 , t>0 

d t p - nAp + (u • V)p = -M 2 u 3 x£K 3 , t>0 

* V • u = iel 3 , t>0 
K u\t=o = u , p\t=o = Po 

where the unknown functions u = u(x,t) = (tti,tt2, u 3 ), p = p(x,t) and p = p(x,t) are the 
fluid velocity, the thermal disturbances and the pressure, respectively. The parameters v > 0, 
n > and g > are the viscosity, the thermal diffusivity and the gravity force, respectively. 
The parameter N > is Brunt- Vaisala frequency (stratification-parameter). Recall that 
A := {31 + d\ + dl)I 3 , V := (81,82, 83) and e 3 := (0, 0, 1). 

Our method follows the ideas based on BMN. For the limit equations, we show that it 
is equivalent to the 2D-Navier-Stokes equation^], which is known to have, in the almost 
periodic setting, a unique global solution, see for example |15j . Then, we show that the 
global existence for the remainder equations in the limit equations. Since we handle not only 
periodic functions, we have to introduce a new analytic functional setting, which is more 
suitable for the almost periodic situation as the second author did for the rotating fluid case 
in [20]. More precisely, a straightforward application of an energy inequality is impossible if 
the initial data is almost periodic. To overcome this difficulty, we use ^ 1 -norm of amplitudes 
with sum closed frequency set. We recall the analytic functional setting (see [20] ) as follows: 

Definition 1.1. (Countable sum closed frequency set.) A countable set A in M 3 is called a 
sum closed frequency set if it satisfies the following properties: 

A = { a + b : a, b G A} and - A = A. 

Remark 1.2. If {ej} 3 =1 is the standard orthogonal basis in M 3 , then the sets Z 3 , {miei + 

\/2m2e2+m 3 e 3 : mi, • • • , m 3 G Z} and {miei+m 2 (ei+e 2 v / 2)+rn3(e 2 + e3v / 3) : mi,m 2 ,m 3 G 
Z} are examples of such countable sum closed frequency sets. Clearly, the case Z 3 corresponds 
to the periodic. Each of the other two cases are dense in M 3 and therefore they correspond 
to "purely" almost periodic setting. 

Definition 1.3. (An i l -type function space) Let BUC be the space of all bounded uniformly 
continuous functions defined in K 3 equipped with the L°°-norm. For a countable sum closed 
frequency set Acl 3 , let 

X A (R 3 ) :={ U =Y^ u n e in ' x G BUC(R 3 ) : u_ n = u* n for n G A, ||u|| := ^ \u n \ < 00}, 

nGA neA 

where n* is the complex conjugate coefficient of u n . 

The second condition in Definition 11.11 is needed to include real-valued almost periodic 
functions in X A . 

Remark 1.4. Note that functions in i 1 do not necessarily decay as x — > 00. Also, this 
almost periodic setting is in general, different from the periodic case since the frequency set 
may have accumulation points. The almost periodic setting is somehow between the periodic 
and the full non-decaying cases. 
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Now, we define anisotropic dilation of the frequency set as follows. 
Definition 1.5. For 7 = (71,72) £ (0,oo) 2 , let 

(1.2) A(7) := {(71^1,72^2,^3) £ M 3 : (ni,n 2 ,n 3 ) £ A}. 

Now, we specify the following Quasi-Geostrophic equation (a part of limiting system), and 
assume (for the moment) that it has a scalar global solution 9 = 6(t) = 9(t,xi,X2,xs), 



(1.3) 



M - A 3 6 + (-A fc )-Va [( w . V ) ((-A fc )V20))] = 0, 

«;= (-a^c-A^-va^a^c-Afc)- 1 ^) 

i 6*(t)| t= o = 6*o, 



where A^ := + <9 2 2 and A3 := d\ t + <9 2 2 + <9 2 3 . We will show that the initial value problem 
for the QG equation admits a global-in-time unique solution in C([0, 00) : X A ) with the 
initial data O = -d X2 {-A h )- l / 2 uo t i + d xl {-^ h )~ x / 2 UQ^. 

More precisely, we give an explicit one-to-one correspondence between the QG and a 2D 
type Navier-Stokes equations (for the existence of the unique global solution to 2D-Navier- 
Stokes equation with almost periodic initial data, see [E]). Now we state our main result. 

Theorem 1.6. Let A be a sum closed frequency set. There exists a set of frequencies dilation 
factors I\A) C (0, oo) 2 such tha$: 

for any 7 £ T, for any zero-mean value and divergence free initial vector field uq £ 
initial thermal disturbance po £ v > 0, k > and T > 0, there exists Nq > g 

depending only on v, k, uq, po such that if \N\ > Nq, then there exists a mild solution to 
the equation (TP]) . u{t) £ C([0,T] : X A ^) with zero-mean value and divergence free, and 
p(t)eC([0,T] :X A ^) . 

Remark 1.7. For the periodic case, we do not need to restrict the frequency set to T i.e. 
we can take T(A) = (0, oo) 2 . However, the computation in this case is more complicated and 
needs a "restricted convolution" type result in the spirit of [2]. 



2. Preliminaries 

Before going any further, we first recall the following few facts about the space X A : 

• (X A , || • ||) is a Banach space, and any almost periodic function u £ X A can be 
decomposed u(x) = T, ne AU n e mx , where each "Fourier coefficient" u n is uniquely 
determined by 

u n = lim t— r / u(x)e lx ' n dx, 

\B\^oo \B\ J B 

and B stands for a ball in M 3 (see for example [?])• 

• X A is closed subspace of FM, the Fourier preimage of the space of all finite Radon 
measures proposed by Giga, Inui, Mahalov and Matsui in 2005 (see |12 |, 113^ [T^]). 

• Leray projection on almost periodic functions P = {Pjk}j,k=i,2,3 is defined as 

Pjk-Sjk + RjRk (l<j,k<2) 
with 5jk is Kronecker's delta and Rj is the Riesz transform defined by 

d_ 

dxj 



R j = —(-A)-2 for j = 1,2,3. 
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the complement set T c is at most countable. 



SLIM IBRAHIM AND TSUYOSHI YONEDA 



The symbol &(Rj) of Rj is irij/\n\, where i = y— T (see [7]). Let P be the extended 
Leray projection with Fourier-multiplier P n = {P n ,ij }i,j=i,2,3,4 given by 



P 



n,ij ■- 



Sij-ffi (1 < *, J < 3), 
<5j,- (otherwise). 



• Helmotz-Leray decomposition is defined on almost periodic functions in the same way 
as in the periodic case. Namely, u is uniquely decomposed as 

U = W + V7T, 

where tt = — (— A)~ 1 div u and w = Pu. 
Now we rewrite the system (jl.ip in a more abstract way. Set N := My/g and v = 
(«i,«2,«3,«4) := (ui,u 2 ,u 3 ,^-p). Then v solves 



(2.1) 



dtv — vAv + NJv + V3P = — (v ■ V 3 )v 
v\t= = v 
V 3 • v = 0, 



with i> = diag(z/, u, v, k), the initial data v = («o,i, «o,2, ^o,3) "a?Po), V 3 := (<9i, <9 2 , <9 3 , 0), 



J 



/'O ^ 



-1 

\o 1 y 



and (w • V3) = (vidi + t> 2 <9 2 + t>3<9 3 ). 

Observe that under the condition N > g we have A/" > y/g and therefore ||f 0,4 1| = || ^Po II < 
\\po\\. We will assume this condition throughout the paper. 

Applying the extended Leray projection P to (|2.ip . we obtain 



(2.2) 



dv/dt + {-DA + NS)v = -P(v ■ V 3 )v, 
v\ t =o = Pv = v Q , 



with 5" := PJP. Recall that for \n\h 7^ 0, the matrix S n := P n JP n has the following Craya- 
Herring orthonormal eigen basis {q^, g" 1 , g°, g„ } (see J3j [9]) associated to the eigenvalues 
{iu n , — iuj n , 0, 0} with 



and 



\n\h 
\n\ 



\ n \h 



n\ + n 2 . 





in 



in 



(,Ql,ni 92,n) l\,ni ol,n) 
,n' 92,n> ^3,n' ^4,n) 

(«l,n 5 92,71) 93,71) ^4 .1(7 

( „div „div div n div\ 
\Hl,ni 92,71) 93,71) 94,7i J 



■(iuj n n 1 n 3 ,iuj n n2n3, -i\n\ 2 h u; n , \n\ 2 h ) 



(-iw n nin 3 , -iiv n n 2 n 3 , i\n\ 2 h uj n , \n\\) 



-1* 



a* 



V^lnl 2 , 
1 -(-n 2 ,m,0,0) = 



PI/7 

1 



n 



(ni,n 2 ,n 3 ,0), 



'Ill 

In 
'in 
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where q^* = (%)* is the conjugate of <jr„. The case when \ri\h = and rc.3 7^ 0, we define 

--(1/2,1/2,0,1/^2) 
=(-1/2, -1/2, 0,1/^2) 
=(-1/a/2,1/V2,0,0) 
=(0,0,1,0). 

In fact, for = and 713 7== 0, we have 5 n = P n JP n == 0. However, the above choice of the 
basis is uniquely determined by the conditions {vq^-q}?) = (^p); (^Qn 1 'In 1 *) = (ht 1 ) an d 
(^<7n'<Zn*) = ^ f° r ^|2,6p . Moreover, the divergence-free condition requires that (v n (t)'Qn V ) = 0> 
giving^ := (0,0,1,0). 

Using Craya-Herring basis, one obtains an explicit representation of the solution to the 
linear version of (|2.2p . For n G A and v n :== (v n ,i, v Uj 2, v n ^, v n> 4) such that v n ■ n == 0, we have 

JNS„ " _ ST^ n croJNS n cr _ ST^ cr ia uj n Nt a 

n ~ / j n e Hn ~ / j u n e Hn 1 

(7 6{-l,0,l} (7 6{-l,0,l} 

with n := (n, 0) = (ni, 712, 713, 0), 

n = ^ <°C and <° :==(£„ •<£<>*)• 
(7oe{-i,o,i} 

Similarly, write a solution u of ()2.2|) as 

neA 

From (|2.2|) . we derive for n G A, 

(2.3) d t v n {t) == -v\n\ 2 v n {t) - S n v n (t) - iP n E m)%(i) with (n • v n (t)) = 0. 

n=fc+m 

In the sequel, we do not distinguish between n and n unless a confusion occurs. For n G A 
we have 

e t7vs ^ n (t)= ^ <°(*)e w5 »e= E < o (*)^ ubW c, 

<7 6{-l,o,i} <r e{-i,o,i} 
where a°°(t) :== (v n (t) ■ q%°*)- From equation (|2.3p . we get for oo = —1,0, 1, 

^<°(i) = -a^(t)((|n| 2 P + iV5 n )C-C*) 

- * E ^C(C' m )M'C). 

n=fc+m, (Ti ,«72S{— 1,0,1} 

Note that P n is self adjoint and P n q%°* = C*- Setting c£°(i) := e- itArCT ° w "<°(i) leads to the 
following equation 

$<#>(<) = -C(t)|n| 2 (i>c-C*) 

- < E e^-c-C(C-^)(C-C*)- 

n=fc+m, cti,CT2S{— 1,0,1} 

where, w^ fcm := (— aou n + o"iWfc + f2W m ). Now we split the nonlinear part into the "resonant" 
(independent of N) and non "resonant" two parts defined by 

B° n °{g°\h°*) :== -i E (C • • C*)CC 2 

and 

B^{Nt t ^,h n ) :== -i E (C -^)(C -C*)5rC 2 exp« fcm iVt), 

' nkm ' 
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respectively. In addition, observe that we have the following estimates: 



| e -H«l 2 *^o(jvt, g ai ,h a2 )\\ < $2\\g ai \\\\h a2 \ 
\e- u \ n \ 2t B°°(g ai ,h a2 )\\ < % I \\g (T1 \\ \\h a2 \\ 



(2.4) 

t l/2 

(for do = —1,0, 1) obtained by estimating the first derivative of the heat kernel as follows 

-v\n\H 



sup 

n6A 



in e 



The constant C v > is independent of N. 

Then we have the following equations: 
(2.5) 

'd t cl{t) = -y\n\ 2 cl{t) + E (v2)ef i,„, lF (B° n (c°\c^) + B° n (Nt,c°\c°z)) 

d t C(t) = - (*±S) |n| 2 <°W + E( CT1 , CT2)6{ -i,o,ip (fl^(c^,c^) + fl-«(iVt, C - 1 , C ^)) , 

for cro = ±1. Prom the condition u^km = 0' we easily see that the terms B^(c , c 1 ), 
^(c- 1 ,^ 1 ), ^(c ,^ 1 ), B^V), B^V), ^^c ,^ 1 ) and B^(c°,c°) disappear. 
Now, we define the "limit equations" by 

'd t cl{t) = -v\n\^l{t) + Bl{c^^) + Bl{c\c- l ) + B^\c l ), 

(2.6) < 

d t C(t) = - (2±*) |n| 2 C(t) + E( CT1 , CT2)e{ _i,o,i }2 \D 5 n (c ,T1 ,c^), a = ±1, 

where L> := {(0, 0), (-1, 0), (0, -1)} for a = 1 and D : = {(0, 0), (1, 0), (0, 1)} for a = -1. 
Formally, we can get (|2,6p from (I2.5P when N ^ oo. We will justify this convergence in 
Lemma 13.21 Now we show that there is more non trivial cancellation in the limit equations. 
More precisely, 

Lemma 2.1. We have 

B°(c\c- 1 ) + B n (c-\c 1 ) = 0. 
Proof. To prove the lemma, it suffices to show 

(2.7) (ql ■ m){q^ ■ q° n *) + {q^ 1 ■ k)(ql ■ q°*) = for any n = k + m with u k = u m . 
First we show that co^ = uj m if and only if 

k,m € {n € 1? : |n|| = An§} for some A > 0. 
(<J=): This direction is clear. Thus we omit it. 

(=^): Rewrite the identity ujk = w m as F(X) = F(Y), where X := \k\1/k^, Y := IH^/m^ 
and F(X) := Xj(X + 1). Since the function F is monotone increasing, we see X = Y. This 
means that 

u j\ k \ h a j_\ m \h 

= vX 7713 = "7^' 

We only consider the case &13 = and 7713 = since the other cases are similar. A 

direct calculation shows that 

(X \ ( -l 0*\ 1 fk h -m h m 3 \k\ \ 

1 -1 u\i x o*n 1 ( k h • m fe fc 3 |m| A . 



V2A|m||fc||n| h , V A Vl + A 2 , 
By fc 3 = ^ and m 3 = Mfc, we have (|277)l . □ 
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Now we show that the function c° in the limit equations satisfies a quasi geostrophic (QG) 
equation type and that this QG equation is equivalent to the 2D type Navier-Stokes equation. 
By the following lemma, we can see that the function c° satisfies the QG equation (|1.3p . 



Lemma 2.2. Let \n\h = \/ n \ + n | • ^ e resonant part B®(c°, c ) can be expressed as follows: 

R°(r° r°) = - V W * rn)\m\ h 

n=k+m 

Proof. Since q Q n := j^(-n 2 , ni, 0, 0) and q® = 1^(1%^ + \Mh<lm) for ra = /c + m, we have 
B° n (c°,c°) = - Y, clcliql-im^ql-qt) 

n=k+m 

= ~ E \hT\ I ( fc2mi ~ kim2 "> ' (I^Ufffc* + l m U«m)) C fc C m 

n=fc+m WM^Ih 

n=fc+ m \ K \h\ n \h 

, , \ n \h 
n=k+m 

Since k x m = — (m x fc), we see that 



£ r^(^m 1 -A; 1 m 2 )(^.^) C fc c^ = 0, 

n=K+m 

which leads to the desired formula. □ 

Now we show that there is a one-to-one correspondence between the QG and a 2D type 
Navier-Stokes equations. 

Lemma 2.3. Let Ah = d 2 1 + <9 2 2 and 



w 



:= (wx(xi,X2,x 3 ,t),w 2 (xi,X2,x 3 ,t)) := £ ^i,n(t)e m ' x , £ w 2 , n {t)e in ' X 

VneA nGA / 



and define 6 = 9(t,xi,x 2 ,xs) := (— Aft) x ' 2 rot 2 w, with rot 2 is the 2 dimensional curl given 
by 

rot 2 w = d 2 w\ — d\w 2 . 
Then, w solves the following 2D type Navier-Stokes equation 



(2.8) 



dtw — Aw + (w ■ V 2 )w + V 2 p = 0, 
V 2 • w = 0, ty|t=o = w 



if and only if 6 solves (|1.3p . 

Proof. Recall that —Aft = —(3^ + <9 2 2 ). First observe that for 9 = (-Aft) _1 / 2 rot 2 u; 
EneA 6n(t)e m ' x , we have by V 2 • w = 0, 

: V ^-§ n (t)e inx , -i V ^-§ n (t)e in - x ) = (d 2 (-A h )9, -d^-A^O). 



i/' 



v n h ^-1 n ft 

nGA nGA 



Then, apply rot 2 to (|2.8p . we get 



(2.9) 



9trot 2 'u; — Arot 2 w + (w • V 2 )rot 2 w = 0, 
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here, we used the fact that 

(2.10) (w ■ V 2 )rot 2 w = rot 2 [(io • V 2 )H- 

Finally, apply (— A^) -1 / 2 to both sides of (|2.9|) . we see that 9 = (—Ah)~ 1 ^ 2 rotw satisfies 
the desired QG equation (|1.3p . Conversely, applying L := (— (— Ah)~ 1 d X2 , (— Af l )~ 1 d xi ) 
(which commutes with A) to (|2.9p . and by f)2.10[) we can see that Lrot 2 is nothing but the 
two dimensional Leray projection. Therefore, this implies (|2.8p as desired. □ 

Remark 2.4. We refer to |15] for the existence of the unique global solution to 2D type 
Navier-Stokes equation (|2.8[) with almost periodic initial data. 

In what follows and in order to show the main theorem, we need the following lemma 
(which is needed only for the almost periodic case) on the dilation of the frequency set 
(|1.2p . This kind of restrictions is technical. However we do not know whether or not such 
constraints are removable. This means that the general almost periodic setting seems to 
remain open. 

Lemma 2.5. For n,k,m G A and 7 = (71,72) G (0,oo) 2 , define n = ('jini, 7 2 n 2 , 77,3), 
^ = (71 72^2,^3) onrfm = (71ml, 72m2, 7713). Let 

PnUl)--=\n\ 8 \k\ 8 H 8 II 

CT e{-i,i} 3 

Given a frequency set A, there is T := T(A) C (0, oo) 2 s.t. for any 7 G T, P n k m {l) 7^ for 
any n, k, m G A such that (rih, kh,mh) 7^ (0,0,0). 

Remark 2.6. If n h , k h , m h = 0, then JB* 1 ^ 1 , c ±x ) = 0. 
Proof. Define T by 

r := {7 G (0,oo) 2 : P n km{l) / for all n,k,m£A with \n\ h , \ k\ h , \m\ h / 0}. 
Nothe that if 7 G T, then cjZ?~ 7^ 0. We show that T cannot be empty. By a direct 

ThnrTYl 

calculation, we have 

Pnkmin) = |n| 8 |fc| 8 |m| 8 

((Wn + ^~ k + 0Jm)(-Un + UJ~ k + Wm)(Wfi - + Wft) (cc% + ^ - U^)) 2 

= \n\ 8 \k\ 8 |m| 8 (a;| - ( WjE + ^) 2 ) 2 " K " ^) 2 ) 2 
= |n| 8 |fe| 8 H 8 ((^- W |-^) 2 -4^) 2 

= |n| 8 |fe| 8 |m| 8 (4 + w| + u4 - 2 W ?ui - 2u&a;£ - 2^?) 2 

I~l4|7|4l~l4 . I ~ |4i 7 14 I ~ 14 I ~ |4| f |4| ~ |4 

= |m| + |n| |re|^|m| + |n| |re| \m\ h 

-2\n\ 2 \h\l\k\ 2 \k\ 2 h \m\ 4 - 2\fL\ 2 \h\ 2 h \k\ 4 \m\ 2 \m\ 2 h - 2\h\ 4 \k\ 2 \k\ 2 h \m\ 2 \m\ 2 h 

o 2 ; 2 26 o 2 j 2 26 

= — 3n 1 K 1 m 1 7 1 — 3n 2 /c 2 m 2 7 2 

o 2r 2 2 4 2 o 2r 2 2 4 2 o 2r 2 2 4 2 

-3n 1 fe 1 m 2 7 1 7 2 - Sn^m^^ - 3n 2 A; 1 m 1 7 1 7 2 

o 2r 2 2 2 4 o 2; 2 2 2 4 o 2; 2 2 2 4, 

-3n 2 fc 2 m 1 7 1 7 2 - 3n 2 fc 1 m 2 7 1 7 2 - 3n 1 fc 2 m 2 7 1 7 2 -\ . 

Since \n\h, \k\h, \m\h / 0, then the highest order terms never disappear. This means that 

I {7 : ^n,k,mPnkm{l) = 0}| = 0. 

Thus the complement set of V is countable (which means that T is a non-empty set). □ 

Now we show that the limit equations have a global solution. In the almost periodic case, 
the non-resonant part B^ 1 (c^ 1 , c^ 1 ) disappears just by restricting the frequencies set to A(7). 
However, the periodic case is more subtle as we need a lemma on restricted convolution (see 
0). 
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Lemma 2.7. Let A be a sum closed frequency set. If A = 1? , take 7 G (0, oo) 2 , otherwise 
we restrict it to 7 € T(A). TTien /or o"o = ±1 £/iere exists a global-in-time unique solution 
c a °(t) to equations l2M such that c a °(t) G C([0,oo) : ^(A(7))) with «°(t) • n) = /or a// 
n G A( 7 ) and cg^i) = 0. 

Proof. Recall that n = (71 ni, 72^2, ^3), k = (71 £7, 72^2, £3) and rh = (71 mi, 72^2, ^13). First 
we consider the almost periodic case. By restricting 7 G T, we can eliminate the worst non- 
linear term using Lemma 12, 51 More precisely, for all n G A and a = (07, o"2, 03) G {— 1,1} 3 , 
the term 

B?°(c°\c°*) = Yl disappears. 

w' T - =0 

fikrh 

Then we have two coupled linear equations for {c~ 1 } n and {c*} n . In this case, the global 
existence will immediately follow from estimates (|2.4p . 

However, the periodic case requires more details. For a > and p > 1, define the weighted 
norm as 

VneA / 

The main step is to show an a-priori bound on c^ 1 in Since the 3D type Navier-Stokes 
equation (|2.6p is subcritical in the space i\ with s > 1/2, then a bootstrap argument (using 
the dissipation, see [17} Proposition 15.1] for example) enables us to conclude that 

c ±1 (t)G^ c ([0,oo):4)nL£((0,oo):^) for 1< 3 < 2, 

whenever the initial data c ±:L (0) G More precisely, by (|2.6p we write the following mild 
formulation: 

(2.11) c ^{t) = e^\ n \ 2t c^{t )- E / e- i ^ ( *- T+to)|n|2 J B^(c <Tl ,c ,T2 )o!T 

K, CT2 )6{-l,0,l} 2 \ J D t( > 

for < to < t. We have a good estimate for the heat kernel with fractional Laplacian, 



\n\ e 



S e -M 2 * = t -s/2\ t l/2 n \s e -\ty2n\z < t -s/2 C e < t _ s / 2 - e /2. ^ 



1 + |tV2 n |iA - 1 + jnpA 



for < t < 1 and e > 0. Then by Holder's and Young's inequality for the discrete case, we 
have the following estimate from (|2.11|) : 



(2.12)|K°(t)|U < C K>L 



K°(io)|U + [ (t-r + to) ' 

Jt 



t 

s/2-e/2 



E (11^^)11,111^(^11,0 + || C -(r)||,o||c-(r 

(<ri,(T2)e{-i,o,i} 2 



If 1 < s < 2, then ||c CT °(f)||^> is finite since the right hand side of (|2.12p is finite. This means 
that c a ° G L^ c ((0,oo) : i^) for 1 < s < 2. Then, thanks to Bernstein's lemma (which can be 
applied only for the periodic case), ^ n \c n \ < \\c\\^ for a > 3/2, we get an a-priori bound 
of the £ 1 -norm. Now we show an a-priori bound of c^ 1 in t\. Multiply equation (j2.6|) by 
l n | 2 ( c n°)* an d summing, we obtain 



E E^(NC(*)) 2 + (\n?C(t)f 



o- e{-i,i} nez» 



EE E K°( C °\c°i)\n\ 2 cr(t). 

nez 3 o- e{-i,i} (o-i, o- 2 )e{-i.o,i} 2 \D 
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Notice that 



(cr 



[{vn ■ qT)e' itNaouJn ]* = (v* ■ q°°)e itNaouJ " 



-ao 



Moreover, from one hand, a direct calculation using Holder's then young's inequalities 
shows that 



(2.13) 



E 5 « 



n 

^'i-'i^n^ + || c ±i|| 4 || c 0||^< e || cT i||| +c 
for sufficiently small e > 0. On the other hand, we show that 



< C\\c Tl \\ e 2 I ||c" ||^|| 



±l112 "^lll + ll^lllllc " 2 



II go IK 



(2.14) 



E^( 



c*,c*)\ n \*<g, 



< 



E E 

n —n=k+m 



l^llcmPllcfc 1 !! 71 !!^ 1 ! 



Using identities a x (a x /?) = — \a\ 2 /3 and (a • (j3 x 7)) = ((a x j3) • 7), we see that 



E E (irn-qf)^ 1 



[in x zn x 



Tl ))) - ±1 - ±1 - Tl 



n —n=k+m 

E (^^fe 1 ) ((9m 1 xiii)'(»xg„ Tl ))cf^c; 

n — n=k+m 

- E E ' ^XSm 1 x ik)(in x ^c^c^cj 1 

n —n=k+m 

- E E ( im • 9fc ^(^m 1 x im)(in x ^^c^c^cj 1 

n —n=k+m 



By the skew-symmetry, namely, 



E E ( im " ^^(^ x im)(in x g* 1 )^ 1 ^ 1 ^ 1 

n —n=k+m 

- E E ( in " 9fc ^(^m 1 x im)(in x gj^c^c^cj 1 

ra —n=k+m 

E E ( in " ^X^m x im)(in x q-J 1 )^ ^cj 1 , 

rn —m=k+n 



the second term disappears which obviously leads to (|2.14p . Now to have the a priori bound 
of c^ 1 in l\, we apply a smoothing via time averaging effect due to [2]. 

Proposition 2.8. Restricted convolution^ Theorem 3.1 and Lemma 3.1] Assume that the 
following holds for \n\h / 



sup 

Tl 

k:k-\-m-\-n=0,k^T l i 



(2.15) 

n 

for every i = 1, 2, • • • , where 

Si := {k :2 i < \k\ < 2 i+1 }, X (n, k, m) 



1 if P nkm (l)=0 
if P nkm (l)^0. 
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Then we have 
(2.16) 



E^( 



-2" "<- 2 

< Cll^lli+ellc 



=1=1 1 1 2 



Now, combining (|2.13|) and (|2.16p . we obtain the following estimate on c in l\, namely. 



(2.17) 



.±i 



| < 11^(0)111 +^^(l|c ±1 («)ll| + l|c ±1 ( S )|||||c°( S )|||) 



ds. 



Moreover we have the following energy inequality: 



V + K 

~~2 



(2.18) 



< 11^(^11,0 + ||c u (t)||,o + 

< 11^(0)115, + ||c°(0)L 2 i . 



"^(s)\U+v\\c°(s)\U)ds 



In fact, multiply the first equation of (|2.6p by c°*(t) and the second one by c ±1 *(t), we 
have (|2.18|) . Since all convection terms disappear due to the skew-symmetry, for example, 



^(cV)- 1 *) = -*E E 



(9fc-m)(^-Q n 1 )4c^c, 1 



n — n=k+m, — w,\+aj^,-)-a;l 1 =0 



E 
E 



E 



n -n=fc+m, -cji +o;i+a;i 1 =0 



E 



m m=k+n, ui n +uh— a) m =0 

= -ElWc- 1 )-^)- 

To apply Gronwall's inequality, we need the following definition. Let C be the positive 
constant appearing in (|2.17|) . Then define h as 

rr+h' ^ 

/i := inf{/i' G [0, oo) : / W^ 1 (s)\\ 2 £l ds < — for any r > 0}. 
./t 2 2G 



Note that h is independent of N and can be chosen positive thaks to (|2.18|) . From (|2.17|) 
and an absorbing argument, we see that 

sup ||c ±1 ( S )||f 1 < 2||c ±1 (0)||f 1 +2C /"* sup ||c ±1 ( S ")|||i||c°( S , )ll|i^ / t < h. 

0<s<t 2 2 JO 0<s"<s> 2 1 



Then by Gronwall's inequality, 



sup W^&Wl <2||c ±1 (0)|L 2 1 expf2G f \\c°(s)\\lds) for t 
•<s<t 2 2 V Jo 1 / 



< h. 



0< 

Iterating the same argument, one more time, we obtain 

sup llc^s)!! 2 ! < 2\\c ±1 (h)\\ 2 fl exp (iC I ||c°(«) ||5ids^ for h<t<2h. 
o<s<t 12 12 V Jh 1 / 

Note that J Q * ||c°(s)|| 2 ids is always finite for any fixed t > 0, since there is a global solution 

to the 2D Navier-Stokes equations in If-type function spaces. Fixing T > and repeating 
this argument finitely many times, we have an d priori bound of c^ 1 in l\ over t G [0,T]. 
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Finally, to use Proposition 12.81 on the restricted convolution, we need to verify that (|2.15p 
holds. Observe that 

P nA . n . k (l) = Ht|fc| 4 |n + A ; | 4 + H 4 |^|n + fc| 4 + H 4 |A ; | 4 |n + ^ 

-2\n\ 2 \n\ 2 h \k\ 2 \k\l\n + k\ 4 - 2\n\ 2 \n\ 2 h \k\ 4 \n + k\ 2 \n + k\\ 

-2\n\ 4 \k\ 2 \k\ 2 h \n + k\ 2 \n + k\ 2 h 

= \n\\kl + l.o.t. 

where l.o.t. stands for lower order terms. Thus, it follows that P n & n _fc(l) is a polynomial 
of degree eight in k% with a nonzero leading coefficient whenever |ni| + \ri2\ 7^ 0. Then for 
fixed k\, &2 and n, there are at most eight ^3 satisfying x( n i k, —n — k) = 1. Thus, 

^2 N^xO,^-™- k) < ^2 Ik^xi^K-n- k) 

2 i <|fc|<2 i + 1 0<|fc| h <2 l + 1 ,fc 3 6K 

^ 8 E E \k\^ <8Y,2 2 ^ +1 h^ <C2\ 

j=l 2J<|fc| h <2i+ 1 j=l 

□ 

3. Proof of the main theorem 

Before proving the main theorem, we first mention the local existence result. Using esti- 
mate (|2.4p . we obtain a local-in-time unique solution to (|2.5p in C([0,T] : ^(A)) as stated in 
the following lemma. 

Lemma 3.1. Assume that c(0) := (0)}neA,<T e{-i,o,i} G ^(A) and c o°(°) = f or °~o G 
{ — 1,0, 1}. Then there is a local-in-time unique solution c(t) £ C([0, Tl] : ^(A)) andcQ°(t) = 
for o"o G { — 1,0, 1} satisfying 

(3.1) r L >— sup ||c(t)|| < 10||c(0)||, 

l|c(o)ir o<*<t l 

where C is a positive constant independent of N. 
Proof. First we recall the mild formulation of ()2.5|) : 

<£(*) = e-l"l 2i c°(0) 

+ ^ / e -"(*-«)M a ^(c CTl ,c ,T2 ) + J B0(^,c ,71 ,c' T2 ))ds 

(<n >0 *)e{-i,o,i} a 

and 

C(t) = e-^l"l 2 *C(0) 

+ J2 I e-^ ( *- s)|n|2 ( J B^(c CT1 ,c ,J2 ) + B^(Nt,c ai ,c U2 )^ ds. 

(eT 1)0 - 2 )e{-i,o,l} 2 

By (|2.4p . we have the estimates 

\\c°n(t)\\ < ||c°(0)|| + C.t 1 / 2 £ f SUP || C <^)II SUP \\d»(8)\\) 

(CT1 , (T2)e{ _ 1 ,0,l }2 V0< S <t 0<s<t 



and 



< IIC(o)|| 

+cw/ /2 Yl f su p ii cCT1 (*)ii su p ii cfT2 ( s )ii) 

, w ,n,12 V0<Kt 0<s<t J 



( CT i l( r2)e{-i,o,i} 
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These d-priori estimates of sup t ||c°(t)|| and sup t ||c CT °(t)|| give us through a standard fixed 
point argument the existence of a local-in-time unique solution (for the detailed computation, 
see [12] for example). □ 

Let b a °(t) be the solution to the limit equations (|2.6p and c a °(t) be the solution to the 
original equation ()2.5|) . The point is to control, in the ^ 1 -norm, the remainder term r°°(t) : = 
c n°(t) ~ ^n°(^) (°o = — 1)0) 1) by the large parameter N. More precisely, r° and r a ° satisfy 

0-2' 



d t r° n (t) = -u\n\ 2 r° n (t) + ]T [B° n (r^ , c^) + B° n {V* , r°>) + B° n (Nt, , d 

K,^2)6{-1,0,1} 2 



and 



W(t) = -(^)|n| 2 C(t) 



(<ri,<T 2 )e{-i,o,i} 2 

respectively. Once we control the remainder term in the ^ 1 -norm, we easily have the main 
result by a usual bootstrapping argument (see [20j for example). Now we show the following 
lemma concerning the smallness of the reminder term. Let b(t) := {bn (t)} n <=A,a a e{-i,o,i} an d 
r(t) := {r£°(i)} n6AiCToe{ _i A i}. 

Lemma 3.2. For all e > 0, there is Nq > such that ||r n (t)|| < e for < t < Tl and 
\N\ > Nq, where Tl is the local existence time (see Lemma \3 . 1\) . 

Proof. To simplify the remainder equation, we introduce the following notation. Let 

K°(r,c,b): = ]T (BZ°(r°\c° 2 ) + B%>(b°i,r° 2 )), 

(<Ti,<T 2 )e{-i,o,i} 2 

K°(Nt,c): = £ B°°(Nt,c°\c^). 

(<ti, CT2 )G{-1,0,1} 2 



We rewrite the remainder equations as follows: 

idtrlit) = -v\n\ 2 r° n (t) + R° n (r,c,b) + R°(Nt,c) 

\d t r^(t) = -(^)\n\ 2 r^(t) + R^(r,c,b) + R^(Nt,c) for a = -l,l. 



(3.2) 



To control r, the key is to control R°(Nt, c) and R°°(Nt, c) in To do so, we need to 

analyze the following oscillatory integral of the non-resonant part as follows: 

n=k+m,u°, ^0 * ^ nfcm 

1 1 rtkm' 

and 

(3.3) H°°(Nt,c) := J] r ^ 1 ,^ 2 ). 

K,^2)6{-1,0,1} 2 

Note that we have the following relation between B and B {R and 7?.): 

9t (Nt, <f 1 X 2 )) = K° (Nt, <f 1 , h°* ) + B%> (Nt, d t g^ , h^ ) 

+B^(Nt,g°\d t h^) 

and 

d t (K^(Nt,c)) = R^(Nt,c)+K^(Nt,d t c). 
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To control r, we split (|3.2f) into two parts: finitely many terms and small (in ^ 1 (A)) 
remainder terms, respectively (cf. [TJ Theorem 6.3]). For 77 = 1, 2, • • • , we choose {sj}'^ =1 C N 
(si < S2 < • • • ) in order to satisfy — V v )r\\ — > (77 — > 00), where 



ni, 



, n St) eAin^ni (k ^ I), \rij\ < 77 for all j = 1, • • • , j. 



The choice of ni • • • n s?; is not uniquely determined, however this does not matter. Then 
we can divide r into two parts: finitely many terms r ni , • • • , r ns ^ and small remainder terms 
{{I ~ V v )r n } neA . 

Remark 3.3. We have the following estimates: 

IIWpy.Tvoiio < ^(i + v 2 ) 1/2 WP v c\\l 

\\T v K a {V v y,cM\ <^ + r} 2 ) 1/2 \\V r) y\\{\\c\\ + \\b\\) 

and 

|||n| 2 ^|| < (l + ^H^yllo 
for < t < Tl {Tl is a local existence time, see (|3.1|) ). where 

P(rj) := maxjlw^l" 1 : k = fa,- ■■k Sv , n = n lr -- ,n Sri , m = n - k}. 

Note that /3(r/) is always finite, since it only have finite combinations for the choice of n, k 
and m. We can also have the same type estimate for ||9f'P r; c|| using (|2,6p . 

We use a change of variables to control R° and R a ° . Let us set y as 

y° n (t) := r° n (t) - K n (Nt,V v c) and y°*(t) := C(i) - K°°(Nt,V v c). 
From (|3.2p . we see that 

d t (y° n + K° n ) = - v \n\\yl + nl) + Rl{y" + n\c,b) 

+R° n (Nt,V v c) + R° n (Nt, (I - V v )c), 

d t (y°°+K ) = -r^^|n| 2 (C+K ) + <°(y CT0 +^,c,6) 

+R a n °(Nt,V v c) + K°(Nt, (I - V v )c). 
Now we control V v y° and V^y 170 for fixed 77. By (|3.3p , 
(3.4) d t V v y° n (t) = -u\n\ 2 V v y n + V v R n (V v y,c,b) + E° n , 

d t V v y%>(t) = -(^^|n| 2 PX°+^n (ny^,6)+^, 

where 

E° n : = -V v iZ n (Nt, V v d t c) + V v R n (V v K° n (Nt, V v c), c, b) 
-u\n\ 2 V v n° n {Nt,V v c) + V v R° n (Nt, (I - V v )c) 



and 



K° = = -V V K° (Nt, V v d t c) + V V K° (V V K°° (Nt, V v c),c, b) 

V + K 



l V v Kl {Nt,Vr,c) + V v R a n °(Nt, (I - Vr,)c) 
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Note that ()3.4p are linear heat type equations with external force E° and E m) . Thus the 
point is to control E° and E a °. By Remark 13.31 we can see that for any e > 0, there is r/o 
and Nq (depending on r/o) such that if N > Nq and r\ > ijq, then H-E^H < e and ||-E' cro || < e. 
Thus we have from (|3.4p . 

\\V v y n (t)\\ < £(^(1 + ^)11^(5)11 

+(1 + ^V'll^WIKIK^II + \\K')\\) + z)ds 

and 

\Wn°(t)\\ < 1^ (c u , K (l + V 2 )\\V ri y^(s)\\ 

+(1 + ^ 2 ) 1/2 H V°(s)\\(\\c(s)\\ + ||6( fl )||) + e)cto. 

By Gronwall's inequality, we have that for any e > 0, there is r]Q and Nq (depending on 
r/o) such that if rj > rjo and N > Nq, then ||"Pr/y°|| < e and \\P v y a ° \\ < e for < t < T^. 
Clearly, we can also control (I — V v )y with sufficiently large r] (independent of N), and 
T , v lZn° (Nt,Vrjc) with sufficiently large N for fixed r]. Thus we can control r for sufficiently 
large r\ and N. 

□ 
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